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ABSTRACT 

We present HST WFPC2 V band imaging for 23 Ultraluminous Infrared Galaxies 
taken from the QDOT redshift survey. The fraction of sources observed to be inter- 
I acting is 87%. Most of the merging systems show a number of compact 'knots', whose 

. colour and brightness differ substantially from their immediate surroundings. Colour 

' maps for nine of the objects show a non-uniform colour structure. Features include 

blue regions located towards the centres of merging systems which are likely to be 
areas of enhanced star formation, and compact red regions which are likely to be 
dust shrouded starbursts or AGN. The host galaxies of the QSOs in the sample were 
Q^i found to be either interacting systems or ellipticals. Our data shows no evidence that 

I ' ULIRGs are a simple transition stage between galaxy mergers and QSOs. We propose 

O I an alternative model for ULIRGs based on the morphologies in our sample and previ- 

ous N-body simulations. Under this model ULIRGs as a class are much more diverse 
^ ' than a simple transition between galaxy merger and QSO. The evolution of IR power 

source and merger morphology in ULIRGs is driven solely by the local environment 
and the morphologies of the merger progenitors. 

' Key words: infrared: galaxies - galaxies: active - galaxies: Seyfert - galaxies: star- 

^ I burst - Quasars: general 



1 INTRODUCTION 

The existence of galaxies that emit more radiation in the 
rest-frame infrared than at all other wavelengths combined 
was esta blished with the first extragalactic mid-IR obser- 
vations ( Low & Kleinmann 1965 ; Rieke fc Low 197^ ) . The 
Infrared Astronomical Satellite (IRAS), launched in 1983, 
was the first observatory with sufficient sensitivity to detect 



above 10^^ in the local {z < 0.5) Universe, outnumber- 
ing Seyferts, Quasi-Stellar Objects (QSOs) and Starburst 
Galaxies of comparable bolometric luminosity. 

A subset of the IRAS galaxy population that has been 
extensively studied are the Ultraluminous Infrared Galax- 
ies (ULIRGs), those with Lir > 10^^ hg^ Lq (or equivalently 
Leo > lO" "/ig5^L0, where Leo is the IRASQOfim flu x ). Re- 



these Lj Urn inouE Infrared Calaxioe in lar go numbore ( Seifa 



et al ig84;-Rn'Bran-Rnbinsnn et al 1984 



Zanders fc Mirabe 



1996). Luminous Infrared Galaxies (LIRGs) are broadly de- 
fined as any object with Li,- > IO^^LqQ These objects be- 
come the dominant extragalactic population at luminosities 



* where Lir spans 8-lOOO^m and is computed using the fiux in all 
four IRAS bands according to the prescription given by Perault 
(1987). 



cent deep submillimetre su rveys (Hughes et al 1995 ; Blain 
et al 1999|; [Lilly et al l"999| [Barger, Cowie fc Sanders 199E| 



Peacock et al 2000) suggest that ULIRGs are an important 



component of the Universe over 2 < z < 7. The evolution 
of ULIRGs, their power source and the trigger mechanism 
behind the IR emission are still the subject of vigorous de- 
bate. Although it is now generally accepted that a mixture 
of star formation and AGN activity powers the IR emission, 
the dominant power source, and how ULIRGs evolve over 
time, are still unanswered questions. There are many simi- 
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Table 1. Ultraluminous QDOT Galaxies Observed By HST 
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Coordinates and V band magnitudes are taken from the HST Planetary Camera images and are in the HST flight filter system. 
^Spectral classification, taken from Lawrence et al 1999 and the NASA/IPAC Extragalactic Database: 'Sb' - Starburst, 'Syl/Sy2' - 
Seyfert 1/2, 'BLRG' - Broad Line Radio Galaxy, 'HX' - High Excitation but not Seyfert 2. ^QOfim luminosity taken from Lawrence et 
al 1999 and converted to Hq = 65, CIq = 1 "^Logarithm of the F606W band luminosity, computed from the F606W band images and in 
units of bolometric solar luminosities (3.823 X 10'^'' ergs s~^). *As measured from the WF3 F814W band images, but see also Borne et 
al 2001, in preparation. ^ IRAS 25/im to 60^m flux ratio. ^First number is from the 7 band system presented in Lawrence et al 198S 
and described in §4.1., second letter is from the 3 band system described in §4.1. ^Asymmetry statistic, computed using Equation ]A 
and described in §4.1 



larities betwe en ULIRGs and starburst galaxies ( Joseph^ 



Wright 11985 
et al 190 



Rowan- Robinson & Crawford 1989; Condon 



Conversely, many ULIRGs displa y nuclear emis- 
sion lines characterist ic of Seyfert galaxies (San ders et a] 
al 199E| ). A study of 15 ULIRGs (|Gcnzel et 



198^ 



I 



UTi et 



al 19981 using the Infrared Space Observatory {ISO) found 
that the IR emission of 75% were predominantly powered by 
star formation and 25% by an AGN, although at least half 
the sources showed evidence for both types of emission. The 
spatial density, bolometric emission and luminosity function 
of ULIRGs are comparable to tho se of Quasars in the local 
Universe. It has been suggested (Sanders et al 1988) that 
ULIRGs represent the dust enshrouded precursors to opti- 
cally selected QSOs and that all QSOs emerge from a lu- 
minous infrared phase. The majority of ULIRGs appear to 
lie in interacting or mergin g systems, although the exact 



fraction is still not certain ( ganders et al 1988 



Leech et al 



1994; [Elements et al 1996). A merger would provide an ideal 
environment for triggering a nuclear starburst or AGN and 
thus would be a possible triggering mechanism for QSOs. A 
merger bet ween two spi rals is also likely to form an ellipti- 
cal galaxy (Barnes 1989), therefore ULIRGs may also have a 



pivotal role to pl ay in elliptical galaxy evolution ( Kormendy 
& Sanc fers 19^ . The relation of ULIRGs to their higher 
luminosity counterparts, Hyperluminous Infrared Galaxies 
(HLIRGs, th ose with Ljr > 1 0^'^ h'^^ Lq) is also not fully 
understood (Hines et al 1995; Rowan- Robinson 200C). Al- 



though HLIRGs as a class appear to be an extrapolation of 
ULIRGs to extreme luminosities, there exists the possibil- 



ity that so me fraction of HLI RGs are an entirely different 
population (Farrah et al 2001). 

In this paper we present Hubble Space Telescope images, 
taken with the F606W filter using the Planetary Camera on 
WFPC2, for 23 ULIRGs and compare them to WFPC2 WF3 
F814W band data for 22 of these sources. Sample selection 
and observations are described in §2. Data reduction, pho- 
tometry, colour map construction and point spread function 
subtraction are described in §3. Results are presented in §4. 
Discussion of the results can be found in §5 and conclusions 
are summarised in §6. 

We have taken Ho = 65 km s~^ Mpc~^, f2o — 1.0 and 
A = unless otherwise stated. 



2 SAMPLE 

2.1 Sample Selection 

The 23 galaxies in our sample are randomly selected from 
those ULIRGs d iscovered in the QDO T all sky /7?/15 galaxy 
redshift survey ( Lawrence et al 1999| ). This survey consists 
of an all sky sample of 2387 IRAS galaxies brighter than 
the IRAS PSC 60/im completeness limit of 0.6Jy, and is 
complete to greater than 98%. All the galaxies in our sample 



have 60/im luminosities in the range 10^^'^^ — W^^'^'^ h^^ Lq. 
The redshift range is 0.14 < 2 < 0.35 with a median redshift 
of z = 0.18 and a mean redshift of z = 0.20. 



© 2001 RAS, MNRAS 000, ??-?? 



HST Imaging of QDOT ULIRGs 






Figure 1. HST F606W images of IRAS 00275-2859 to IRAS 06361-6217. 
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Figure 2. HST F606W images of IRAS 06561+1902 to IRAS 14337-4134. 
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HST Imaging of QDOT ULIRGs 





Figure 3. HST F606W images of IRAS 16159-0402 to IRAS 20109-3003. 
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Figure 4. HST F606W images of IRAS 20176-4756 to IRAS 23220+2919. 
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2.2 Observations 

The data were taken in cycle 8 between October 1998 and 
July 1999 using the Wide Field Planetary Camera 2. The 
coordinates of each source were centred on the Planetary 
Camera CCD, selected for its superior pixel scale (0.046" 
pixel"^ as opposed to ~ 0.1" pixel"^ for the WFC) and full 
sampling of the PSF over the Wide Field Camera CCD's. At 
the mean redshift of our sample the pixel scale of the Plane- 
tary Camera corresponds to a physical size of 0.2kpc. Expo- 
sures were taken using the F606W filter, which corresponds 
approximately to the V band filter in the Cousins system. 
The observations were performed in 'Snapshot' mode, each 
object being imaged for 600 seconds. Each observation was 
split into 2 exposures of 300s to facilitate the subtraction of 
cosmic rays. 

For 22/23 of the sample, data were taken using the Wide 
Field Camera 3 CCD and the wide band F814W fiher, which 
corresponds closely to the / band filter in the Cousins sys- 
tem, as part of a separate program (prop. ID 6346, PI K. 
Borne). These observations were also taken in 'Snapshot' 
mode, each object being imaged for 800s. 

The conversion between HST filter system magnitudes 
and Cousins magnitudes is not straightforward. Depending 
on the filter used for the observations the applied conversion 
factor can be strongly related to the object spectrum. For 
the F606W filter the dependence of the conversion factor to 
Cousins V band on the object spectrum is quite large, and 
can introduce photometric errors of Am = 0.2 or more. For 
the conversion between F814W magnitudes and Cousins / 
band the dependence is much smaller, generally Am < 0.05. 
In order to avoid these associated uncertainties, magnitudes 
are given in the H ST filter system using the Vegamag syn- 
thetic zeropoints (Holtzman et al 1995), and have not been 



converted into the Cousins filter system. The differences be- 
tween our F814W magnitudes and the equivalent Cousins / 
band magnitudes are likely to be very small, with Am < 0.1. 
For the F606W filter, which is both wider and redder than 
the V band Cousins filter, the match is much poorer. Hence 
we do not make comparisons between our F606W band mag- 
nitudes and previously quoted V band magnitudes, but do 
make comparisons between our F814W band magnitudes 
and previously quoted Cousins / band magnitudes. 



3 DATA ANALYSIS 
3.1 Data Reduction 

Upon arrival all datasets were recalibrated with the best 
available reference files, using the IRAF task CALWP2, and 
combined into a single image using the IRAF task Crrej. 
Following sky subtraction, warm pixels were removed by lin- 
ear interpolation from their immediate neighbours. Photo- 
metric solutions were calculated using the Synphot pack- 
age from the STScI, rece ntly u pdated to match the solutions 
given by Holtzman et al (1995). For those QSOs in the sam- 



ple that were saturated in the central regions magnitudes 
and fluxes were computed using PSF fitting photometry. 
Corrections were made for detector gain (for these obser- 
vations the gain was 7e~ /DN) and aperture size. Absolute 
magnitudes were derived using the expression: 

-1/2. 



Figure 5. Histogram showing the number of sources per band 
in the 7 band classification system of Lawrence et al 1989. A 
summary of this classification scheme can be found in §4.1. 



We estimate that our derived relative magnitudes have an 
associated photometric error of approximately 4%. Fluxes 
in Jy were obtained by multiplying the fiux in ergs cm~^ 
s~^ by the central wavelength of the F814W filter in 
angstroms, and then dividing by the central wavelength in 
Hz. As such they do not include any bandpass correction. 

No k corrections have been applied to the calculated ab- 
solute magnitudes. Trial k corrections for both the F606W 
and F814W filters wer e computed using the spectral syn the- 
sis package Pegase (Fioc & Rocca-Volmerange 1997). At 



M 



25 - 5log[{2c/W0h)(l + z){l 



1 + ^ 



(1) 



the mean redshift of the sample, z — 0.2, the k correction in 
the F606W band was found to be of the order Am = 0.15 
for a lO" year old spiral galaxy, with the correction in the 
F814W filter being slightly lower. Even at the highest red- 
shift in the sample, z = 0.345, the k correction was not much 
greater than this (Am = 0.18). It was found that, if the age 
of the system was varied across 10^ — W^" years and/or a 
different source spectrum was assumed, then the computed 
k correction varied across the range —0.02 < Am< 0.3. The 
errors introduced in applying these synthetic k corrections 
therefore probably significantly outweigh any corrective ef- 
fect. There are also no applied corrections for reddening due 
to dust extinction. 



3.2 Colour Maps 

Colour maps are an invaluable tool for highlighting mor- 
phological features not apparent in single band images. Con- 
structing colour maps for the F606W and F814W band data 
presented here is not straightforward as the data were taken 
on different CCD's. The Point Spread Function (PSF) for 
any of the WFPC2 CCD's is strongly dependent on posi- 
tion, wavelength, and CCD characteristics. These differences 
must be carefully accounted for when matching the two im- 
ages. 

The F606W and F814W band images were first scaled 
to a common exposure time and gain. Although the readout 
noise does not scale linearly for the PC and WF3 chips, this 



© 2001 RAS, MNRAS 000, ??-?? 



8 D. Farrah et al. 



effect is not significant when the exposure time of one im- 
age is at least 10% of the other. Following the advice of the 
STScI, the images were then deconvolved wi th appropr iate 
PSFs using the Richardson-Lucy algorithm (Lucy 1974) as 
implemen ted in the IRAF task LuCY. The Tiny TIM v5.0 
software ( Krist 1994 ) was used to generate PSFs for the ap- 
propriate filter and chip positions. The four frames in the 
F606W and F814W observations were then combined into 
single images using the IRAF task Wmosaic, which cor- 
rects for geometric distortion, offsets and scale differences 
between the Planetary Camera and the Wide Field Cam- 
eras. This task also degrades the Planetary Camera V band 
images to the resolution of the Wide Field Cameras, effec- 
tively convolving the PC image. Any residual PSF mismatch 
is thus due to differences in the effective PSF's of the images 
achieved by the initial deconvolutions. Since this PSF mis- 
match could in principle be significant we performed sim- 
ulations to assess the accuracy of this procedure by con- 
structing colour maps of stars visible in both the PC and 
WFC CCD's. It was found that the maps of these stars had 
a featureless colour structure. Although some level of PSF 
mismatch is undoubtedly still present it appears to be negli- 
gible compared to noise in the original images or large angle 
scattering within the WFPC2 optics. 

A further potential source of extraneous colour struc- 
ture with this method is that of deconvolution artifacts. 
These can include sharp colour peaks caused by amplifica- 
tion of poisson noise in the original images, or very bright 
point sources surrounded by dark rings which are caused by 
large numbers of iterations on a bright point source within 
an envelope of extended emission. To minimise the occur- 
rence of these artifacts we examined the images after each 
iteration. It was found that deconvolution artifacts only be- 
came visible several iterations after convergence between the 
original image and the deconvolved image had been achieved 
(measured via a reduced statistic). As we halted the de- 
convolution immediately upon attaining convergence, the ef- 
fect of artifacts on the achieved colour structure is insignif- 
icant. 



3.3 QSO Host Galaxies 

Host galaxies can only be resolved in QSOs if light from 
the central bright source, and from the host galaxy, can be 
separated. The most convenient way to accomplish this is 
by fitting a point source template to the QSO; the template 
being established from a set of observations of stellar PSFs 
that are a good match in colour, filter and detector position 
to the targets. For the optical QSOs in this sample no addi- 
tional observations were made of nearby stars, and no suit- 
able observed PSFs were available in the HST PSF Archive. 
Synthetic PSFs were therefore generated using Tiny TIM. 

Ten times oversampled PSFs were generated. These 
PSFs were then shifted by non-integer pixel distances, re- 
binned to Planetary Camera resolution and convolved with 
the PC pixel scattering function. The agreement between 
these synthetic PSFs and observed (stellar) PSFs is gener- 
ally excellent within radii of about 2". Beyond this radius 
scattered light in the WFPC2 optics can have a significant 
effect, an effect that is not modelled by the Tiny TIM soft- 
ware. The effect of any systematic uncertainties in the PSF 



have also not been examined, such effects are however likely 
to be negligible. 

The resulting PSFs were then centroided and normal- 
ized to the image via a reduced fit- It was found that 
all the QSO images were saturated in the central regions 
to varying extents, mostly confined to the central 3x3 pix- 
els and never beyond the central 5x5 pixels. In all cases the 
centralmost possible pixels were used, avoiding those that 
were saturated or that lay in the diffraction spikes. This 
method allows the image and PSF to be registered to within 
~0.2 pixels (estimated by visual inspection), and also gives 
a starting value for the PSF normalisation. 



3.4 Profile Fitting 

In order to determine the morphology of the host galaxies of 
the QSOs in the sample, both de Vaucoleurs and exponential 
disk profiles were fitted to the surface brightness profiles of 
the host galaxies. The de Vaucoleurs profile is of the form: 



/(r) 



J -7.67[(r/re)l/*-l] 



(2) 



and the exponential disk profile is of the form: 

/(r) = Joe-''-'"'' (3) 

One-dimensional surface brightness profiles of the 
sources were extracted by fitting annuli spaced at 1-pixel 
intervals, fixing the ellipticity (at e — 0.05), semimajor axis 
and position angle for each fit. The PSF profile together with 
either a radial de Vaucoleur or disk profile were then fitted 
to the source surface brightness profile, treating the PSF 
normalisation as a third free parameter. This helps to min- 
imise as far as possible the systematic errors introduced by 
subtracting off a 'best guess' PSF and fitting a galaxy pro- 
file to the remaining light distribution. Fitting was carried 
out using Levenberg-Marquardt least-squares minimization 
in order to robustly derive the best fit parameters. Two di- 
mensional isophote fitting, using the IRAF task ELLIPSE, 
was also used to investigate the ellipticities and isophotal 
twisting (large changes in azimuthal angle as a function of 
radial position) for each host by iteratively fitting isophotes 
of constant surface brightness whilst allowing the ellipticity 
and position angle to vary. 

In order to quantify whether a host galaxy was detected 
or not we imposed the condition that a PSF -I- galaxy to 
QSO fit must be better than a pure PSF fit to at least 
99% confidence, computed via an F test. We excluded pix- 
els within the central 0.2" as the PSF within this radius is 
highly undersampled. We also excluded all pixels beyond 2" 
because of light scattered within the WFPC2 optics. 



4 RESULTS 

4.1 Morphology and Colours 

Positions, redshifts, magnitudes, spectral types, luminosities 
and merger classification for the sample can be found in Ta- 
ble ^. The HST Planetary Camera images of the sample can 
be found in Figures |l[ § and ^ For the two binary sys- 
tems where the galaxies are clearly separated the magnitude 
of both sources together is given. Magnitudes of individual 
galaxies in binary systems can be found in Table ^ 
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Figure 6. Comparison between the Lawrence et al classifica- 
tion sciieme, tiie qualitative classification scheme described in §4.1 
and the asymmetry statistic A described in §4.1. The three dif- 
ferent systems, although based on different criteria, overall give 
similar descriptions. The three systems are therefore measuring 
related physical quantities, as expected. The circles denote the 
mean asymmetry for each Lawrence et al band. 



Contributions to the measured magnitudes from emis- 
sion lines could in principle be significant, as many ULIRGs 
are kno wn to be strong emission line galaxies ( Veilleux, Kim 
& Sand 3rs 199!: ). At the mean redshift of our sample the Ha 
lies near the centre of the F814W bandpass, and the OIII 
line lies nearly in the centre of the F606W bandpass. We 
evaluated uncertainties in our measured magnitudes by com- 
parison with Hot and OIII luminosities taken from Veilleux, 
Kim & Sanders ( |l999| ). The F606W band luminosities of our 
sample span the range 10.0 < Ly < 11.4, with all but three 



having Lv > 10.3 (where Lv is the logarithm of the F606W 
band luminosity in units of bolometric solar luminosities). 
The F814W band luminosities, although not quoted in Table 
|l[ span a similar range. In contrast, the dereddened Ha lu- 
minosities from Veilleux, Kim & Sanders (199!;) lie in the 
range 7.0 < Lna < 10.5, with most lying in the range 
8.0 < Lhc < 9.3. Similarly, the OIII luminosities span 
the range 7.0 < Loin < 10.0, with most falling between 
7.0 < Loin < 8.0. The uncertainty in our quoted magni- 
tudes and colours due to these emission lines is therefore 
almost certainly much less than Am= 0.1 for all the objects 
in the sample. 

In order to provide a quantitative description of the 
morphologies in the sample we have classified the sources 
according to the seven band prescription given by Lawrence 
et al (1989): 

0: Source has no companions within 200Kpc and shows no 
signs of interaction or merger. 

1: Source has a faint companion (between 4 and 2 magni- 
tudes fainter than the source) between 40Kpc and 200Kpc 
away. 

2: Source has a bright companion (less than 2 magnitudes 
fainter than the source) between 40Kpc and 200Kpc away. 
3: Source has a faint companion less than 40Kpc away but 
shows no signs of interaction. 



4: Source has a bright companion less than 40Kpc away but 
shows no signs of interaction. 

5: Source is interacting with companion, and there are signs 
of tails, loops, or bridges. 

6: Source is merging; either there is obvious disturbance 
and /or there are two or more nuclei in a common envelope. 

No source in the sample is an isolated undisturbed spiral 
or elliptical galaxy. Using the above classification scheme, 
3/23 of the sample are classified as not interacting (groups 
0-2), 3/23 of the sample are classified as having a nearby 
companion but showing no signs of interaction (groups 3 - 
4), and 17/23 are classified as interacting or peculiar (groups 
5-6). Figure ^ shows a histogram of the number of sources 
in each band. If it is assumed that the three sources with a 
nearby companion are interacting then this gives a fraction 
of ULIRGs observed to be interacting of 87%. If it is assumed 
that the sources with a nearby companion are not interacting 
then this gives a minimum observed interacting fraction of 
74%. 

For some purposes it is informative to classify the 
ULIRGs into broader categories based solely on qualitative 
descriptions of their optical morphology. We have therefore 
classified our sample as follows: 

Class A: There are two separate sources, but close 
enough to reasonably assume that they are interacting, or 
they are interacting but still physically distinguishable as 
separate objects. 

Class B: The source is not distinguishable into separate 
objects and there is no QSO activity, but there are clear 
signs of interaction such as tails, bridges and and multiple 
nuclei. 

Class C: The source contains an optical QSO with or 
without signs of ongoing interaction. 

These classes are chosen such that sources can be classi- 
fied unambiguously. They are not intended to imply an evo- 
lutionary sequence. By these criteria, 5 sources are identified 
as class A, 14 as class B and 4 as class C. The classification 
of individual objects is given in Table 0. Analysis using this 
qualitative classification scheme is given in §5. 

A quantified estimate of the degree of disturbance in an 
object can be obtained by computing some form of asym- 
metry statistic. Numerous different methods exist, a popular 
approach is to compute the 180° statistic A. Here we use a 
modified form of the equation presented in Br inchmann et 
al ( |199^ ) and used by e.g. Serjeant et al ( |200C| ): 



A = 



2V.. 



(4) 



where lij is the original image, 7^^ is the image following a 
180° rotation, and k is the sky asymmetry statistic: 



(5) 



where Gij and Gij are areas of sky equal in area to the ob- 
ject image. There are two difi'erences between the formula 
used by Brinchmann et al and our formula. The extra factor 
of 2 in both equations ensures that the asymmetry statistic 
is normalised to between (perfect 180° symmetry) and 1 
(perfect 180° antisymmetry). The modulus on the denomi- 
nator ensures that sky values are not overestimated. For the 
QSOs in the sample the asymmetries are computed for the 
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Table 2. Binary System magnitudes 



Source 


RA 


Dec 


My 


M; 


02587-4336 A 


59 43.07 


25 


05.43 


-20.50 


21.50 


02587-4336 B 


59 43.48 


25 


00.88 


-21.62 


22.81 


20176-4756 A 


21 10.67 


47 


08.12 


-20.97 


22.16 


20176-4756 B 


21 11.44 


47 


09.16 


-21.45 


22.52 


Coordinates are 


J2000 and 


are 


given 


as mm ss.ss 





host galaxies rather than for the QSOs themselves as the 
QSOs will otherwise dominate and lead to unrealistically 
low values. For each object the asymmetry was calculated 
for all possible positions within the object and the lowest 
value was then selected. The computed asymmetry values 
can be found in Table ^. 

A comparison between the Lawrence et al classifica- 
tion scheme, the qualitative classification scheme described 
above, and the computed asymmetries can be found in Fig- 
ure ^ It is apparent that, despite the different criteria used 
in these three systems to classify the objects in the sample, 
the asymmetry statistic A gives similar results to both of the 
more qualitative classification schemes. The three different 
statistical measures are therefore measuring related physical 



quantit es, as expected. 

With the enhanced resolution of Hti'l' it is intorma- 
tive to examine the morphology of the sample in more 
detail. Two of the sample (IRAS 02054-f 0835 and IRAS 
10026-1-4347) are QSOs that show no apparent signs of in- 
teraction. Two objects (IRAS 00275-2859 and IRAS 20037- 



1547) are mergers containing QSO activity. Three sources 
(IRAS 06268+35U9, IRAS 18520-5048 & IRAS 20253-3757) 
are merging galaxy pairs. One source (IRAS 06361-6217) 
resembles a coUisional ring galaxy. 13/23 sources show ev- 



idence :or tidal tails, plumes, and other signs of ongoing 



interaction. In addition these sources all show at least one, 
and in most cases several, compact bright regions or 'knots'. 
Two sources (IRAS 02587-6336 & IRAS 20176-4756) are bi- 
nary systems. Discussion of the morphology of individual 
sources can be found in §4.4. 

Three sources from our sample can be identified as hav- 
ing 'warm' infrared colours (/25//60 > 0.2), whereas ten 
sources have 'cool' infrared colours (/25//60 < 0.2). The re- 
maining ten sources have IRAS upper limits that do not 
allow us to distinguish between these two possibilities. Of 
the three 'warm' sources, two are QSOs, one of which shows 
signs of interaction. Of the ten 'cool' sources, one is an inter- 
acting system containing QSO activity, one is a close pair, 
and the rest are interacting systems. Two of these interact- 
ing systems (IRAS 13469-^5833 & IRAS 18580-^6527) show 
evidence for a violent merger, with large numbers of bright 
knots and large tidal tails. 

4.2 Colour Maps 

Colour maps have been constructed for those nine sources 
that show sufficient spatial extension in both the F606W 
and F814W band images, and are presented in Figure ^ 
Sky background regions, and those regions with low signal 
to noise, have been masked out. The slightly lower S/N in 
the F606W band means that excessively reddened (dark) 
regions with low S/N are masked off. This effect is most 
apparent in IRAS 185804-6527 and IRAS 20109-3003, where 



the maps become increasingly dark towards the edges of the 
objects before being masked off. 

One object (IRAS 00275-2859) displays a smooth colour 
structure (discounting the slightly aberrated centre), per- 
haps indicating the presence of a uniform dust shroud or 
screen obscuring the central illuminating source. The other 
colour maps all display evidence for a large scale non- 
uniform colour structure which cannot be attributed to ar- 
tifacts of amplified noise spikes in the deconvolved images. 
Two objects (IRAS 04384-4848 & IRAS 14337-4134) show 
a central, extended, diffuse blue region. One object (IRAS 
06268+3509) shows several red, compact regions coincid- 
ing with the centres of the two merging galaxies. We inter- 
pret these as heavily obscured star formation regions. The 
remaining sources show a number of blue or red compact 
'knots', in most cases located towards the central regions. 
The blue regions are most plausibly interpreted as areas of 
enhanced, unobscured star formation, and the red regions 
as either dust shrouded star formation or remnant nuclei 
from the merger progenitors. Care must be taken in inter- 
preting such regions from imaging data alone. An alterna- 
tive interpretation is that these regions are reflec ted light 
from an o bscured AGN, as observed in Mrk 463 (Uomoto 
et al 1993| ). Although most ULIRGs (> 80%) are starburst 
dominated, AGN reflected light could still be a potential 
source of uncertainty in classifying knots as most starburst 
dominated ULIRGs have a signiflcant AGN component, es- 
pecially at higher luminosities. Another possi bility could be 



shock fronts from a n expanding superwind (Heckman, Ar- 
mus & Miley 1990). Although a superwind might be ex- 



pected to produce a systematic colour gradient as opposed 
to a point-like stru cture comparisons between loca l ly ob- 
served sup erwinds (Heckman, Armus & Miley 1990; Wang 
et al 1997] ) show that, at the mean redshift of our sample. 



a superwind could create a small structure only a few pix- 
els (~ 5) across. A flnal possibility is that of small holes 
in locally obscured regions, caused by photodestruction of 
dust grains by hot young stars. These alternate possibilities 
are impossible to completely discount without spectroscopic 
data. 



4.3 Host Galaxies 

Four of the objects in the sample contain some form of op- 
tical QSO. Host galaxy images for these four QSOs are pre- 
sented in Figure pi Data for the host galaxies are presented 
in Table |[ 

Host galaxies were clearly detected in all four QSOs. 
For IRAS 20037-1547 & IRAS 00275-2859 both spiral and 
elliptical light profiles were very poor fits to the underlying 
light distribution and the host galaxies are clearly interact- 
ing. For IRAS 02054-^0835 & IRAS 10026+4347, the host 
galaxies were well fitted by an elliptical profile (x^ ~ 1-4), 
with a spiral profile being a much poorer flt to the data 
ix^ ~ 4.0). We investigated the possibility that we were 
merely fltting the bulges of spiral host galaxies by fltting si- 
multaneously elliptical and spiral components in the central 
and outer regions respectively, the radius of each compo- 
nent being allowed to vary. In both cases it was found that 
a 'pure' elliptical profile was strongly preferred. 

IRAS 02054+0835 has a host galaxy which displays 
no evidence for ongoing interactions. The ellipticity is low 
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Figure 7. _F606iy — F81AW colour maps for nine of the 22 sources. The colourbar gives the range in colour for each source. The 
orientation of each map is the same as the source image. The F606W band data has been degraded to F814W band resolution. The 
central region of IRAS 00275-2859 is aberrated due to the very bright central point source. 



(e ~ 0.05) and there is no isophotal twisting. The other re- 
solved host galaxy (IRAS 10026-1-4347), ahhough well fitted 
by an elliptical profile, also shows moderate disturbance. Al- 
though having only two resolved hosts means comparisons 
with other samples are necessarily limited, our host galaxy 
properties are simila r to Early type galaxies in the Forna x 
and Virgo clusters ( Gaon, Capaccioli fc D'Onofrio 1994 ), 
and to Elliptical galaxies at higher redshifts ( |Kelson et a] 
200C|)~1 

Of the remaining 19 objects in the sample without an 
optical QSO, 17 are too disturbed to attempt fitting any 
form of galaxy profile. The remaining two objects (IRAS 
17431-5157 & IRAS 23140-1-0348) have a sufliciently regu- 
lar morphology such that it might be expected that fitting 
galaxy profiles would yield informative results. Although one 
profile was clearly preferred in each case it was not possi- 
ble to obtain a 'good' fit for any profile. Reduced val- 
ues were of the order 3 for the best fit in each case. We 



therefore conclude that these galaxies, despite being regular 
in appearance, have alight profile not consistent with a dy- 
namically relaxed spiral or elliptical galaxy and are therefore 
disturbed. 



4.4 Individual Sources 

In the following descriptions, comments regarding HST im- 
ages refer to the F606W Planetary Camera images. 

IRAS 00275-2859: This IR luminous QSO was discov- 
ered by Vader and Simon (1987). Clements et al <\l996) 



calculate a dust temperature of 54K and a dust mass of 
7.3 X IO^^Mq for this object, assuming a linear emissivity - 
frequency relation and a monolithic dust population whose 
temperature is solely defined by the 60/100/im fiux ratio. 
Zhenlong et al (1991) classify this object as having a near 
(< 10") faint companion but not interacting. Clements et 



al (1996) classify this object as 'disturbed'. The HST im- 
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age shows strong interactions and a double nucleus. The 
northeastern nucleus contains a QSO. There is also a spa- 
tially (but not luminally) symmetric halo extending about 
2". There is one very bright 'knot' near the QSO, and a num- 
ber of dimmer knots near the southwestern nucleus. Canalizo 



Table 3. QSO Host Galaxy Properties 



& Stockton (2001) classify this second nucleus as a giant HII 
region rather than an AGN. The faint companion described 
by Zhenlong et al is revealed to be an undisturbed spiral. 

IRAS 02054+0835: A QSO with no signs of ongoing in- 
teraction. There are, however, five small, faint sources within 
10". 



IRAS 02587-6336: Previous optical imaging ( jZhenlong 
et al 19^) classified this system as a close pair with interac- 
tion. The HST image shows a binary system. Both sources 
show signs of disturbance, although there are no tails or 
loops and the sources appear still physically separate. The 
brightest source is a disturbed spiral, with a smaller source 
to the southwest. 

IRAS 04384-4848: A small, disturbed source with no 
QSO activity. There are a number of compact bright regions 
and two faint tails. 

IRAS 06268+3509: A pair of merging spirals. The 
northern spiral contains several compact bright 'knots'. 
There are four further small sources within 12". 



IRAS 06361-6217: Zhenlong et al ( |l99l| ) classify this 
object as being non-interacting, but with a nearby compan- 
ion. The HST image shows a source with a collisional ring 
galaxy morphology. The south eastern nucleus is very bright, 
with several nearby knots. A faint, elliptical halo surrounds 
these regions, with a semimajor axis length of about 8". The 
western region is smaller and dimmer, with no bright knots. 
There is a small, faint companion 10" to the north, and an 
object 15" to the south which is either a QSO or foreground 
star. 

IRAS 06561+1902: An interacting pair with a close 
companion. Although still physically distinct, the two in- 
teracting galaxies are connected by a faint bridge of nebu- 
losity. This source is striking in that there is a similar group 
of three objects 10" to the southwest, one of which is a QSO. 

IRAS 07381+3215: A smaU, pecuhar galaxy with some 
bright starforming regions. Notable are the two foreground 
stars 15" south east of the source. 

IRAS 10026+4347: This object shows less than 1% po- 
larisation in the optical/UV, rapid X-ray variability and no 
optical absorption, indicative of a low quantity of dusty gas 
in the object and zero intrinsic polarisation, and probably 
indicative of a d irect, unobscured line of sight to an AGN 
disk or nucleus (Grupe et al 1998). The optical and X-ray 



pro perti es of this object are discussed extensively by Xia et 
al (1999), who describe the source as a post-merger system 



with strong Fell emission. It shows a steep X-ray contin- 
uum slope, an X-ray luminosity comparable to that of Syl's 
or QSO's, large rapid X-ray variability, and optical contin- 
uum variability. The intense Fell emission is cited as being 
the result of gas enrichment by supernovae. The i/ST image 
shows a QSO with no apparent signs of interaction. There 
are three very faint sources within 12". 

IRAS 10579+0438: The HST image reveals a small 
source with two bright regions separated by about 1.5", 
linked by a small tail. The western region is brighter, larger 
and more disturbed than the eastern region. The only object 
within 30" is a foreground star. 



Source 


My 


re (Kpc) 


le 


Notes 


00275-2859 








Interacting 


02054-f0835 


-22.63 


3.7 


20.7 


Elliptical, low e 


10026+4347 


-21.14 


6.6 


22.2 


Elliptical, disturbed 


20037-1547 








Interacting 



Properties of the host galaxies of the 4 objects in the sample 
containing an optical QSO. Data for the two interacting hosts 
are not given. 



IRAS 13469 +5833: Ground based optical imaging 
(Leech et al 1994) showed an interacting system with giant 



tails. The HST image shows a large and complex morphol- 
ogy. There is a large, bright, circular region to the west with 
a radius of approximately 1". To the east of this there is 
a cluster of smaller 'knotlike' regions. There are two large 
tails on the east and west sides of the source. Large numbers 
of bright 'knots' are also apparent throughout the source. 
There are two faint sources 15" to the west. 



IRAS 14337-4134: Zhenlong et al (|l99l|) identify this 
source as a close pair with a tidal tail, and classify it as 
interacting. The HST image shows a merging system with 
several bright regions and two small tidal tails. There are 
four bright sources within 25". 

IRAS 16159-0402: This interacting source shows a 
bright, peanut shaped region with some small knots imme- 
diately to the north. Further north there is a larger, dimmer 
region. A tail or plume stretches south to an inclined disklike 
structure. 

IRAS 17431-5157: A bright, disturbed core surrounded 
by a dim asymetric halo. There are a large number of other 
sources within the PC field of view. 



IRAS 18520-5048: Zhenlong et al ( [l99l[ ) describe this 
object as being non-interacti ng bu t having a close faint com- 
panion. Clements & Baker (1996) however argue that this 
system in fact possesses two nuclei separated by 2.5". The 
-ffSr image reveals an interacting pair of spirals. The two spi- 
rals are similar in size and brightness, both having a number 
of small, bright 'knots'. 

IRAS 18580+6527: T his source is cited as a possible 
multip le nucleus ULIRG (Auriere et al 1996). Leech et al 
(1994) tentatively identify this source as a merger system 
interacting with a companion. The western source has a Sy2 
spectrum and the eastern source has an HII spectrum. The 
-ffST image shows a source that is strongly interacting. There 
is a large, bright region to the west of the source, with several 
smaller bright regions surrounding it. To the east there are 
two smaller bright regions, together with a very large tail 
extending approximately 10" to the northwest. The western 
region appears brighter and more active than the eastern 
region. 

IRAS 20037-1547: An interacting system containing a 
QSO. Remnants of one of the merger progenitors lie to the 
east. Physical merging appears to have only just started in 
this system. 

IRAS 20109-3003: A pair of merging spirals. There are 
two compact, bright regions separated by about 1", which 
can plausibly be identified as the galaxy nuclei. The north- 
western region is both larger and brighter than its south- 
western sounterpart. A long 'tail' extends about 10" to the 
west, and there are no companion sources. 
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Figure 8. HST F606W images of the ULIRG QSO host galaxies. Residual PSF features have been masked out. 



IRAS 20176-4756: This source is difficult to interpret. 
The most plausible picture is a binary system consisting of 
a very compact double nucleus source and a slightly dis- 
turbed spiral with a small companion, separated by 30kpc. 
Under this interpretation, this source would be in the earliest 
stages of interaction. It is however not possible to draw solid 
conclusions without redshifts for both sources. The quoted 
redshift is for the spiral. 

IRAS 20253-3757: An interacting pair with no nearby 
companions. The northern galaxy is a disturbed spiral, with 
arms extending to the north and south. A number of lumi- 
nous 'knots' are apparent in the southern galaxy, extending 
east-west, 



the nucleus and the other restricted to elongated or 'fan- 
like' areas. This implies that significant amounts of radia- 
tion are escaping at large angles to the radio axis. The core 
is described as having a compact steep radio spectrum with 
extensions to the east and west. The optical spectrum of the 
associated galaxy is characteristic of young A type stars. An 
intensity contour map of the source shows circular isophotes, 
whereas polarised intensity isophotes are very anisotropic, 
showing a 'fan' to the southwest and a slight enhancement 
to the northeast. The radio axis of the source does not coin- 
cide with any of the polarisation axes of the source. The flux 
f rom the radio core of this object is uncertain. Early st udies 



(lUlvestad 1985) do not ag ree with later VL A stud ies ( Har 



IRAS 20507-5412: This sourcti is small, luminous and 
displays a number of compact bright regions. There are two 
bright sources within 15", to the south east is a spiral and 
to the north east is a small QSO. 

IRAS 23140-1-0348: Yee and Oke (1978) reported broad 
wings on both Ha and H/3 for this radio-loud source. An 
optical polarisation study (Draper, Scarrott & Tadhuntei 
iggaini^veais two polarisation components, one intrinsic to 



vanek fc Hardcastle 1998 ). Heckman et al (1994) identify 



this source as a narrow line radio galaxy. Interestingly, this 
source is a member of one of the 'compact groups' identified 
by Hickson (1982), these groups being characterised by pos- 
sessing at least three members, having a high density and 
a low velocity dispersion. Members of such groups would 
therefore undergo almost continuous gravitational pertur- 
bation and could be candidate sites for triggering ULIRG 
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activity (Borne et al 2000). The ifST image shows a bright 



source with small, dim knots to the east and north. There 
is a dim asymmetric halo, extending towards the northwest. 

IRAS 23220+2919: A small source with a bright nucleus 
to the southeast and an extension to the southwest which 
probably contains a second smaller, dimmer nucleus. There 
are no companion sources within a 15" radius. 



5 DISCUSSION 

5.1 Sample Morphology 

The large scale morphology of ULIRGs has been studied 
by m any a uthors, often with conflicting results. Sanders 
et al ( |l98^ ) found that 100% of a sample of 10 U LIRGs 
showed signs of interaction. Lawrence et al (1989) found 
that, although 19/41 IRAS galaxies with Leo > W^^Lq 
were i nteracting fwhere Lm is the GOfim luminosity com- 



puted i :om the IRAS 60/jm flux), only 2/6 IRAS galaxies 



with Leo > W^^Lq were i nteracting. Near I R imaging of 



a sample of nine ULIRGs (Carico et al 199C) showed that 
four of the sample possessed double nuclei with high activity 
levels, and they estimated the mean starburst lifetime to be 
4 X 10* ye ars. VLA observations of two co mplete samples of 
ULIRGs (sopp, Alexander & Riley 199C) found several bi- 



nary starbursts, and that starburst activity can be triggered 
in both, eit her or neither galaxy in an interacting pair. Zhen- 
long et al (1991) derived an interacting fraction of 61% for 



their sample of 41 ULIRGs. Leech et al (1994) imaged 35 
ULIRGs and found that 23/35 showed signs of interaction. 
Their results al so rev ised the interaction fraction quo ted by 
Lawrence et al ( |1989[ ) from 2/6 to 4/6. Clements et al ( |l99e| ) 
derive an interacting fraction of 90% for their sample of 60 
ULIRGs. They explain the discrepancy between their re- 
sults and previous work by saying that deep optical imaging 
is needed to detect faint signs of interactions and that previ- 
ous stu dies mav have missed such sig n^. Qbservatioii s of the 



IRAS l |jv sample of ULIRGs (iKim fc Sanders 199^) foJ^ 
a luminosity function of the form^^(L)~or^~'^ , strong 
density evolution with a — 7.6 ±3.2, and no evidence for 
clustering. Murphy et al (1996), using optical and near-IR 
imaging of a sample of 56 ULIRGs, quote 95% as inv olved 
in current or recent interactions. Surace et al (1998) per- 
formed deep B and / band imaging using HST for a sample 
of nine 'warm' ULIRGs and found that all showed signs of 
ongoing interactions. They also found that most of their 
sample contained a number of compact blue 'knots' of star 
formation that do not contribute significantly to the bolo- 
metric emission, and that a small number of these knots 
were plausibly active nuclei. Overall they found that their 
results were consistent with 'warm' ULIRGs representing 
a critical transition phase between a galaxy merger and a 
QSO. A complementary study of a sample of 'cool' ULIRGs 
( Surace, Sanders fc Evans 200C| ) found that all of the sample 
were either advanced mergers, or early stage mergers with 
evidence for separate (> 600pc) nuclei. Emission from the 
central regions was, with one exception, cons istent with high 
rates of star formation. Murphy et al (2001), using integral 



field spectroscopy of a sample of 4 ULIRGs, propose that the 
ULIRG phase may be bimodal in time with some ULIRGs 
being early stage mergers and some being late stage merg- 
ers. They further propose that multiple starburst events may 



take place during the ULIRG phase and that some LIRGs 
may be 'resting' ULIRGs. 

The number of sources observed to be interacting in 
our sample of ULIRGs is 87%, with a lower bound of 74%. 
These results are in excellent agreement with the more re- 
cent ULIRG morphology studies using deep, high resolution 
ground and space based imaging, although care must be 
taken in comparing this fraction to lower resolution sam- 
ples. Five of the sources are binary systems where the pro- 
genitors have yet to coalesce. The merging systems show a 
wide range of morphological features and asymmetries, most 
containing at least one compact luminous 'knot'. Of the four 
QSOs two (IRAS 02054-f-0835 & IRAS 10026+4347) show 
no signs of recent interaction. The other two QSOs are still 
interacting. In one case (IRAS 00275-2859) the merger is ad- 
vanced, whereas in the other (IRAS 20037-1547) the QSO 
and galaxy have only just started to physically merge. 

The Sanders et al picture of ULIRG evolution ( ganders 
et al 198^ ) postulates that ULIRGs are the dust shrouded 



precursors to optically selected QSOs. According to this pic- 
ture (hereafter referred to as the S88 picture), interactions 
and mergers between gas rich spirals transport gas to the 
central regions of the galaxies. This large central gas concen- 
tration triggers starburst activity, and in the latter stages of 
the merger creates or commences the fuelling of a central ac- 
tive supermassive black hole. This central black hole rapidly 
comes to dominate the IR luminosity of the ULIRG. In the 
last stages the dust screen shrouding the black hole is blown 
away and the ULIRG evolves into an optical QSO. 

The sample presented here provides excellent evidence 
that ULIRGs are strongly linked with both galaxy interac- 
tions and (to a lesser extent) QSO activity. That QSO activ- 
ity is physically connected with ULIRGs is clearly demon- 
strated by the number of optical QSOs in the sample. If the 
two phenomena were unconnected then the number of ob- 
jects in the comoving volume out to z = 0.4 which are by 



chance b oth ULIRGs and QSOs is 0.1 (Canalizo & Stock- 



ton 2001). The sample includes still physically separate but 



interacting spirals, violent mergers with large amounts of as- 
sociated star formation, QSOs with signs of ongoing interac- 
tion, and QSOs with no apparent signs of interaction. Those 
sources that are classified as isolated and non-interacting ac- 
cording to the system of Lawrence et al are all QSOs. 

Looking solely at the optical morphologies of the sample 
it is, however, far from clear that the S88 picture describes 
the sample well. It is not possible to rule out from this data 
alone that, for the two QSOs with no signs of ongoing in- 
teraction, the IR emission was not triggered by interactions, 
a cou nterpart of which exists for the HLIRG population as 
well ( Farrah et al 2001 ). In two sources, a QSO is coeval 
with galaxy interactions. If it is assumed that the e-folding 
timescale of an Eddington luminosity black hole with accre- 
tion ef ficiency e — 0.1 is a good estimate for the lifetime of a 
QSO (Rees 1984) then it is entirely plausible that for these 
two sources the QSO will burn out before the merger is com- 
pleted and the system is dynamically relaxed. A straightfor- 
ward interpretation of all ULIRGs being a transition stage 
between galaxy mergers and QSOs is not supported by this 
data. 

Previous studies have also highlighted potential prob- 
lems wit h the S88 picture. M id-IR spectroscopy of 62 
ULIRGs (Rigopoulou et al 199£) found that the majority 
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Figure 9. Luminosity in the V band, Ly, versus asymmetry 
factor (after subtraction of QSO in case of sources shown in Fig|8|). 
The different symbols correspond to the classification of section 
4.1 (class A - two separate sources; class B - interacting galaxies; 
class C - QSO). 



were powered by starbursts, and that most did not contain 
a buried QSO. There was no evidence that advanced merg- 
ers are more AGN-like, based on the ULIRG nuclear sep- 
arations. That AGN activity becomes more prevalent with 
merger advancement is a natural prediction of the S88 pic- 



ture. A mid-IR spectroscopic survey of 15 ULIRGs (Genzel 
et al ISps]) found that ^ 80% were powered predominantly 
by starbursts and that ~ 20% were powered by AGN, al- 
though at least half the sample showed clear evidence for 
both types of activity. There was no detected trend for 
the AGN-hke systems to reside in the more compact (and 
hence more advanced merger) systems. Canalizo & Stockton 
(2001), in a study of a sample of 9 QSOs with FIR colours 
intermediate between ULIRGs and QSOs, proposed a model 
in which some ULIRGs evolve to become classical QSOs in 
a short time (< SOOMyr) and that some QSOs are born un- 
der the same conditions as ULIRGs but with a very short 
(< SOOMyr) lifetime. They further propose that the QSO 
nuclear regions are initially shrouded in a dust cocoon, an 
idea qualitatively similar to the S88 picture. 

Our data does not allow us to discriminate between IR 
emission due to starbursts and that due to AGN activity. 
We can, however, examine this topic in alternative ways. If 
the IR luminosities of the ULIRGs with a QSO are com- 
pared to the IR luminosities of the ULIRGs without a QSO 
via a Kolmogorov-Smirnov (K-S) test, the the probability P 
that the IR luminosities of these two populations are drawn 
from the same distribution is found to be P = 0.0411. Thus 
it seems very likely that there is a difference in the distri- 
butions of the IR emission between those sources with a 



in more advanced mergers, a result in direct contradiction 
of the S88 picture and indirectly supportive of Rigopoulou 
et al (1999). A plot of Leo vs. A can be found in Figure ^| 
Amongst those sources with no QSO there is no correlation 
between asymmetry and IR luminosity. Both class A and 
class B objects are found across a wide range of IR lumi- 
nosities with no apparent trends. There is some support for 
the idea that ULIR Gs become more AGN-like at the high- 
est IR luminosities (Veilleux et al 1995; Shier, Rieke & Rieke 



1996 ) as three of the four most IR luminous objects in the 



sample are QSOs. 

We performed simulations to assess the impact of both 
large scale features such as tidal tails, and small scale fea- 
tures such as starforming knots, on our derived asymmetry 
values. It was found that, although both large and small 
scale features contributed to the derived asymmetry val- 
ues, the asymmetries were determined predominantly by the 
large scale morphologies. The asymmetry values are there- 
fore an accurate measure of the degree of disturbance in 
our sample, rather than solely a measure of the number of 
starforming knots or other small scale structures. 

Unfortunately there are no other galaxy samples for 
which direct comparisons between asymmetry values can be 
made. The closest samples available are t he optical star- 
burst sample of Brinchmann et al (1998) and the sub- 
ultralu minous radio selected starbursts of Serjeant et al 
( [2000| ), both of which lie in a similar redshift range, and were 
imaged with the Planetary Camera on HST but using the 
F814W filter. Although the difference in filters means that 
different stellar populations are being sampled, the closeness 
in wavelength between the F606W and F814W filters means 
that this difference should not significantly affect large scale 
morphologies. After accounting for the extra factor of two in 
our version of the asymmetry statistic, our computed value s 
are very similar to those computed by Serjeant et al (2000), 
with a slightly higher mean value. They are much highe r 
than those values presented by Brinchmann et al (1998). 
This implies increasing asymmetry with increasing star for- 
mation rate. 

The images presented here have demonstrated the effi- 
cacy of HST for detecting signs of interaction where lower 
resolution ground based images have either shown no signs of 
interaction or have shown ambiguous results. Three sources 
(IRAS 00275-2859, IRAS 06361-6217 & IRAS 18520-5048) 
have, using ground based imaging, been classified as non- 
interacting or ambiguous. The very high resolution HST 
images of these sources have however shown clear signs of 
interaction, implying that ground based surveys of ULIRGs 
may underestimate the interacting fraction by 15% or more. 

Studies of samples of 'warm' and 'cool' ULIRGs have, 
quite recently, become active ar eas of research. ffff T imaging 
of a sample of warm ULIRGs (Burace et al. 1996) revealed 



that all the sources lay in advanced merger system s. Imag- 



ing of a complementary sample of cool ULIRGs (Surace 



QSO a id those without a QSO (as also found by numer- 



ous previous studies), consistent with some form of evolu- 
tion between a starburst and AGN phase. A plot of Ly vs. 
asymmetry statistic A (Equation ^ can be found in Figure 
^ The sources that contain a QSO have comparable asym- 
metries to the non-QSOs in the sample, with (A) = 0.36 
for the QSOs as opposed to (A) — 0.35 for the non-QSOs. 
This suggests that QSOs (and hence AGN) are not found 



Sanders & Evans 2000 ) revealed that all of the sample were 



either advanced mergers or 'premergers' where the galaxies 
were still separate. It has been suggested that warm ULIRGs 
may represent an important transition stage between cool 
ULIRGs and QSOs. Two of the three sources in our sample 
with warm colours are QSOs, consistent with warm ULIRGs 
containing a higher fraction of QSOs than cool ULIRGs. 
The cool sample however contains one QSO with ongoing 
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Figure 10. 60 micron luminosity plotted against Asymmetry 
for all 23 ULIRGs. Description of the classes A, B and C can be 
found in §4.1 



minous radio quiet quasars (Bahcall, Kirhakos & Schneider 



199(:) show that both quasars reside in apparently normal 
host galaxi es, one being an elliptical and one a spiral. Disney 
et al (1995) present //ST imaging for a sample of four QSOs, 
finding that all four have elliptical hosts. Although they find 
that the morphologies of the hosts are all featureless, they 
argue that the presence of multiple very close companions 
is evidence for interactions as the trigger for the quasar ac- 
tivity. 

We resolved host galaxy morphologies in only two of 
the four QSOs in the sample, finding one to be a smooth 
elliptical and one to be an elliptical with signs of distur- 
bance. The hosts of the other two sources were detected at 
> 99% confidence but it proved impossible to fit an ellip- 
tical or spiral profile to the underlying light distributions. 
The PSF subtracted images show these hosts to be inter- 
acting. Although discussion based on only four sources can 
only be very limited, these results do not contradict the S88 
picture, and are also supportive of a link between merging 
spirals and emerging ellipticals. 



interactions. The mean asymmetry for the warm ULIRGs 
in our sample is 0.404, whereas the mean asymmetry for the 
cool ULIRGs is lower, at 0.376. Although our sample sup- 
ports QSO activity being more prevalent in warm ULIRGs 
than in cool ones, our data does not support the hypothesis 
that warm ULIRGs are more advanced mergers than cool 
ULIRGs. 

Comparison between the colour maps and the F606W 
band images shows that, in most cases, bright regions in 
the F606W band images correspond with very blue or very 
red regions in the colour maps. For three objects however 
(IRAS 13469-^5833, IRAS 04384-4848 & IRAS 14337-4134) 
this is not the case. These objects show blue regions that do 
not coincide with bright regions in the F606W band images 
or the F814W band images. As these regions are located 
towards the centres of the sources, this result can be inter- 
preted as a generally slightly higher star formation rate, a 
local underdensity of dust, or a local IMF skewed towards 
high mass stars in what are probably the most turbulent and 
gas rich regions of ULIRGs. These optically bluer regions 
may also therefore correspond to sites of buried starbursts 
and/or AGN. 



5.2 Host Galaxies 

Previous surveys of QSOs using HST have generally been 
extremely successful in detecting host galaxies. Urry et al 
(1999) give results for a sample of BL Lacertae objects, de- 
tecting host galaxies to a redshift of 2; = 0.7. They find 
that the hosts are all extremely luminous ellipticals with 
{Mi) = —24.6. The derived ellipticities are generally low 
and morphologies are smooth. Imaging of a sample of 14 
quasars (6 radio loud, 5 radio quiet, 3 radi o qui et ultra- 
luminous infrared quasars) by Boyce et al (1995) showed 



that all the radio loud quasars had elliptical hosts, the radio 
quiet quasars possessed both elliptical and spiral hosts, and 
that the radio quiet ultraluminous infrared quasars lay in 
violently interacting systems. Host galaxies were on average 
0.8 mag brighter than L* . Results from imaging of two lu- 



5.3 Optical Colours 

A colour magnitude diagram (Meoe — Msu vs. M814) is pre- 
sented in Figure |l^. The average F814W band magnitude of 
those sources which do not contain a QSO is Mg,iA = —22.9. 
This is 0.3 magnitudes less bright than the field ellipti- 
cal galaxy magnitude in the Cousins / band, M'j — —23.2 



([Pozzetti Bruzual fc Zamorani 1996 ; Efstathiou, Ellis & Pe- 



terson 1988). Still physically separate but interacting sys- 



tems (class A) have a similar optical magnitude to those 
systems that are physically merging (class B), with the class 
A objects having an average magnitude of Msi4 = —22.93 
and the class B objects having Mgu = —22.86. The QSOs 
(class C) are much brighter than the field galaxy magni- 
tude, with Af8i4 ~ —24.4, and are plotted for comparison 
only. Similarly, the two objects without an optical QSO but 
with broad lines (18580-1-6527 & 23140-1-0348) are plotted 
only for comparison. 

Evolutionary tracks for three star formation scenarios 
have been plotted to estimate the contribution of stellar pop- 
ulations of various ages to the optical emission from the class 
A/B obj ects. These tracks were computed using the Pegase 
package (Fioc & Rocca-Volmerange 1997), suitably modified 
to produce output in the HST fiight filter system. An initial 
mass normalization was selected by assuming that the Milky 
Way is a suitable template. The total mass of the Milky 
Way is current ly thought to be about IO^^Mq (Wilkinson 



&i Evans 1999). A rough estimate for the mass of stars for 



two merger progenitors (excluding any dark matter compo- 
nent) is therefore about 2 x 10^^ Mq. Dust extinction effects 
are included in the models for a spherically symmetric sys- 
tem and are computed individually for each timestep. The 
calculated optical depth varies between < r < 35. 

The Av = 1 dereddening vector, plotted on the lower 
left of Figure ^ lies almost parallel to the three evolution- 
ary tracks. Coupled with the uncertainty in the mass nor- 
malisation, this means that only upper limits can be set on 
the ages of the stellar populations. Similarly it is not possible 
to draw conclusions about which star formation scenario is 
applicable, or to construct more precise constraints on the 
ages of the stellar populations, as data is only available in 
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two bands which are quite close in wavelength. These limi- 
tations considered, the optical emission from nearly all the 
class A/B ULIRGs is best described by stellar populations 
of several Gyr or less in age. The most likely interpretation 
is that the optical emission is dominated by old stellar popu- 
lations from the merger progenitors , a result consistent with 
previous studies ( Surace et al. 1995 ) , although a heavily dust 



shrouded starburst cannot be ruled out. 



5.4 Pair Mergers vs. Multiple Mergers 

The number of objects involved in a typical ultralumirious 
galaxy merger can have a fundamental impact on both the 
star formation rate and (to a lesser extent) the merger mor- 
phology. Although the physical mechanisms behind star- 
burst and AGN activity are likely to be no different in pair 
mergers and in multiple mergers, the fraction of ULIRGs 
in multiple merger systems can give a good estimate of the 
fract ion of ULIRGs in compact groups of galaxies (Bekki 
2001), as mergers between more than two ga laxies are ex- 



pected to occ ur naturally in compact groups (Barnes 1989 



Hickson 1997). N umerical sim ulations of a merger between 
five disk galaxies (Bekki 2001) showed that repetitive mas- 
sive starbursts are characteristic of multiple galaxy mergers, 
in contrast to pair mergers where only one or two starbursts 
occur. Mult iple galaxy merging has been previously linked 
to Arp 22 (Diamond et al 1989). M odels for multiple galaxy 
mergers ( TaniguchTTT'Shioya 199^ ) propose that the exis- 
tence of more than two AGN or starburst nuclei is evidence 



for a mult iple merger history. It has been suggested (Borne 
et al 200C) that Hickson Compact Groups, where multiple 
mergers may be expected to take place, may be the progen- 



itors of ULIRGs. It has also been suggested ( Murphy et a] 



2001) that multiple starbursts, as may be expected to occur 



in mergers between more than two galaxies, may be common 
in ULIRGs. 

For those objects in our sample where the progenitor 
galaxies are no longer distinct (class B/C) it is difficult to 
draw reliable conclusions about the number of merger pro- 
genitors, as the complex geometries of gas and dust in these 
systems can easily give rise to multiple bright regions. Even 
using optical spectroscopy to identify separate AGN regions 
as a method for counting the progenitors could be ambigu- 
ous. AGN regions could be heavily obscured by dust making 
them invisible in the optical. Alternativ ely some models sug- 
gest (Taniguchi, Ikeuchi & Shioya 1999) that an AGN could 



form in a galaxy merger via the collapse of a star cluster con- 
taining compact supernova remnants, without the need for 
a 'seed' black hole from the merger progenitors. Although 
we cannot use our data to establish accurate fractions of 
ULIRGs in pair mergers and multiple mergers, we can exam- 
ine whether pair merging and multiple merging does in fact 
take place. For all the ULIRGs in class A the merger appears 
to be taking place between only two galaxies. This is most 
apparent in IRAS 06268+3509 and IRAS 18520-5048. Con- 
versely, some ULIRGs in the sample (e.g. IRAS 13469-1-5833 
and IRAS 18580+6527) show more than two large, well sep- 
arated bright regions which cannot plausibly be interpreted 
as having arisen due to dust obscurati on patterns. Accord- 
ing t o multiple galaxy merger models (Taniguchi & Shioya 
1995 ) this can only arise in mergers between more than two 



can be both pair mergers and multiple mergers, suggesting 
that ULIRGs can be formed both in compact groups and in 
isolated pairs. 



5.5 ULIRG Evolution 

5.5.1 Merger, Starburst, & ULIRG lifetimes 

A prima facie picture of how long each stage in a merger is 
ultraluminous in the IR can be gained from looking at the 
number of sources in the sample in each merger stage. In 
five of the sources the galaxies are still physically separate, 
or they are interacting but still physically distinguishable as 
separate objects (class A). In 14 sources (class B) separate 
galaxies are not distinguishable and there are clear signs of 
ongoing interaction, such as tidal tails, loops, and multiple 
nuclei. In the remaining 4 sources an optical QSO is visible. 
It is therefore plausiblethat physical mergers are the most 
prevalent ULIRG phase, although the optical QSO phase 
and the IR luminous stage before the galaxies are physically 
merging in ULIRGs are both significant. 

Using N-body simulations to assess galaxy-galaxy in- 
t eractions and mergers ha s proven to be an invaluable tool 
(Barnes & Hernquist 1992). The course of a merger event is 



affected by many variables, including (but not limited to), 
angle of approach , relative velocity, relative d isk inclinations, 
and gas masses (Barnes & Hernquist 1996). Hence a com- 



plete discussion will not be presented here. We instead refer 
the reader to Barnes & Hernquist (1992). As an example 
we describ e the N-body simulations of Dubinski, Mihos & 
Hernquist ( 199£| ), specifically their models A and D. Model 
A considers a direct coplanar collision between two equal 
mass galaxies with very low mass dark matter haloes. After 
the first close aproach at t = the galaxies reach a max- 
imum separation of ~ 50kpc at t — 2.5 x 10** years. After 
the second close approach at t 6 x 10* years the merger 
proceeds rapidly and is completed by f ~ 9 x 10* years. The 
same model but with high mass dark matter haloes (model 
D) yields a similar pattern of evolution but with a longer 
merger time of t ~ 2 x 10^ years. 

Overall, N-body simulations of galaxy mergers predict 
merger times (from first close approach to coalescence) of 
approximately 10^ years. The lifetime of a starburst even t 



has been estimated by several authors. Thornley et al (200C) 
find that starburst events are relatively short lived, with a 
lifetime range of 10® — lO'^ years, and that galactic super- 
winds produced by supernova e ma y be responsible for the 
short duration. Genzel et al (1998) derive a starburst age 
range of 10^ - 10* years. Evolutio nary models ap plied to 
interacting and merging galaxies (Bernloehr 1993) give a 



starburst lifetime of 2 x 10^ years, but with gaps between 
starburst events of up to 2 x 10* years. The lifetime of a 
ULIRG can be estimated indirectly in two ways, both inde- 
pendent of the power source behind the IR emission. The 
space density of ULIRGs is approximately the same as that 
of optically selected QSOs at the same bolome tric luminos- 
ity in the local Universe (Sanders et al 1989). If it is as- 



progenitor galaxies. From this, we conclude that ULIRGs 



sumed that both ULIRGs and QSOs are an event in a host 
galaxy largely independent of environment then the range 
in ULIRG lifetimes should be approximately the same as 
that of QSO lifetimes. E stimating the lifetimes of Q SOs is 
difficult. Recent results (Martini fc Weinberg 200C) give a 
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Figure 11. Colour-Magnitude diagram for the QDOT ULIRGs. Evolutionary tracks are for a Salpeter IMF and run from lOMyr to 
lOGyr. The dereddening vector is shown in the bottom left. Time direction for the tracks is L-R. 'Burst' - instantaneous burst of star 
formation at t = 0. 'Exponential' - star formation law of the form SFR = e(-*/iOOOVlOOO. 'Gas Mass' - SFR = Mgas/3000. The class 
C points are not QSO subtracted, and are plotted for comparison only. 



lifetime range of 10''-10** years. Alternatively the maximum 
ULIRG lifetime can be estimated by comparing the mor- 
phologies in our sample to the N-body simulations discussed 
above. The maximum observed separation in the sample is 
30kpc. Assuming that these galaxies are receding from each 



other a 'ter their first close approach gives a time remaining 
to merger of 1.6 x 10^ years (model D) or 8.1 x 10** years 
(model A). Overall therefore we derive a lifetime range for a 
ULIRG of approximately 1 x 10* — 1 x lO^years. This range 
is in broad agreement with previous estimates (Murphy et 



al 1996|j, which quote a range of 2 x 10* — 2 x 10 years but 
argue that the real hfetime probably lies towards the lower 
end of this range. 

Clearly there is some discrepancy between these life- 
times. It seems unlikely that ultraluminous IR activity can 
persist for the entire duration of a merger event, and yet our 
sample show ULIRGs at all stages of a galaxy merger. Sim- 
ilarly it seems implausible that a single starburst event can 
supply the necessary IR luminosity for the entire lifetime of 
all ULIRGs. On the other hand an AGN cannot be responsi- 



ble for the IR emission for the majority of the ULIRG phase 
as most ULIRGs are powered predominan tly by starbursts 
( [Genzel et al 199i ; Rigopoulou et al 199!;). To achieve con- 
sistency multiple starburst events, similar to those pre dicted 
by nu merical simulations of multi ple galaxy mergers (Bekki 



2001) or by previous observations (Murphy et al 2001), must 
therefore be invoked before the final merger. Alternatively 
the evolution of ULIRGs may be more complex than that 
proposed by Sanders et al, or follow a different pattern. 



5.5.2 Starburst & QSO Triggering 

Starburst triggering in galaxy mergers has previously been 
examined using numerical simulations where star formation 



is assumed to follow a density dependent Schm idt law (Mr 
hos & Hernquist 1994; Mihos & Hernquist 1996). In a merger 



between two bulgeless disk/halo galaxies it was found that 
at first perigalacticon bars were formed in each galaxy due 
to tidal forces, causing a rapid inflow of gas into the central 
regions. This central gas concentration triggered powerful 
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starbursts as the galaxies reached their maximum separa- 
tion. This starburst lasts for approximately 10* years, dying 
out well before the galaxies finally merge at ~ 10® years. 
The ensuing starburst during the merger is very weak. In 
a merger between two disk galaxies with a significant bulge 
component the bulges stabilize the galaxies against gas in- 
flow, hence no bars are formed. The resulting starbursts trig- 
gered whilst the galaxies are still separate are therefore much 
less intense than those in bulgeless disks. When the galaxies 
merge however very large quantities of gas are driven to the 
central regions of the merger, triggering an extremely in- 
tense starburst lasting ~ 5 x lO'' years. A few Gyr after this 
starburst has ended, the merger remnant strongly resembles 
a peculiar elliptical galaxy. 

Remarkably, these results depend less on orbital dy- 
namics, such as approach angle and relative disk orienta- 
tion, than progenitor galaxy morphology. The role of or- 
bital dynamics in triggering starbursts is unclear, although 
the relative orientations of the galaxy disks and the peri- 
centric separation after the first close approach appear im- 
portant. Direct encounters produce rapid inflows whereas 
highly inclined passages do not, implying that strong ax- 
isymm etric forc es are needed for prompt inflows (Barnes & 
Hernquist 1996). The predicted star formation rates also im- 



ply that an AGN is not needed to explain the IR emission 
in ULIRGs. Based on numerical simulations and the then 
availa ble optical imaging of ULIRGs, Mihos & Hernquist 
(1996) concluded that most ULIRGs are mergers between 
disk/bulge/halo systems, with the few ULIRGs composed 
of widely separated galaxies being mergers between gas rich 
disk/halo galaxies. 

Observations su ggest that QSOs contain black holes 
with mass > 10* A/0 (Rees 1984). Spiral galaxies are however 



thought to contain black holes with mass < IO^Mq (|Salucci| 
et al 20p0h. Although the masses of black h oles in ellipticals 



cover the range lO^M© < Mbh < 10^°Mq ([Magorrian et a]| 



199^he quantities of gas and dust in such systems are too 
small to trigger ultraluminous IR activity in a merger. Hence 
a > 10* Mq black hole must be formed in a galaxy merger if 
ULIRGs are the transition stage between mergers and QSOs. 
The formation of quasar nuclei in ULI RGs has been investi- 
gated by Taniguchi Ikeuchi & Shioya (1999). They conclude 
that a supermassive black hole (SMBH) of mass > 10*Mq 
can form if either one of the progenitor galaxies contains a 
'seed' SMBH of mass > IO^A/q undergoing efficient Bondi 
type gas accretion over ~ 10* years. 



5.5.3 A New Picture of ULIRG Evolution 

It is apparent from this and previous work that there ex- 
ist serious problems with the S88 picture. Previous results 
have highlighted the ambiguity of whether AGN fraction in- 
creases with advancing merger stage. Our results are not 
consistent with ULIRGs being a straightforward transition 
stage between galaxy mergers and QSOs. Although we find 
2/3 warm ULIRGs contain QSOs, our results are not con- 
sistent with warm ULIRGs being a transition stage between 
cool ULIRGs and QSOs. 

Our results in conjunction wit h previous IR spec- 
trosco py and N-body simulations ([Barnes fc Hernquist 



I992I; iMihos fc Hernnuist 1994 [Mihos fc Hernquist 199^ 



Rigopoulou et al 19991) lead us to suggest an alternative to 



the S88 picture. We propose that ULIRGs as a class do not 
represent a simple transition stage between galaxy mergers 
and QSOs, and that there exists no sequence from 'cool' 
ULIRG to 'warm' ULIRG to QSO. Instead we propose that 
ULIRGs are a much more diverse class of object where the 
time evolution in IR power source and merger morphology 
is driven solely by the morphologies of the merger progeni- 
tors and the local environment. As has been noted by pre- 
vious authors the presence and size of a bulge component in 
merger progenitors has a fundamental effect on the timing 
of gas inflow to the galaxy nucleus. Galaxies with small or 
non-existent bulge components would starburst early in the 
merger, galaxies with a large bulge component would star- 
burst late in the merger. It is apparent from our imaging 
that ultraluminous activity is not mostly confined to late 
stage mergers and that there exist significant numbers of 
ULIRGs at all stages in galaxy mergers, including ~ 20% in 
systems where the galaxies are still separate. This implies 
that there is a substantial variation in the bulge component 
over all ULIRG progenitors. This morphological diversity 
amongst the merger progenitors will give rise to a continuum 
of starburst and AGN power sources as a function of merger 
stage across the ULIRG population. The plausible evidence 
presented here for diversity in environment and hence num- 
ber of merger progenitors (e.g. pair mergers between field 
galaxies and multiple mergers between galaxies in compact 
groups) means that this continuum is more complex than 
may be expected from just pair mergers. It suggests that 
there may be a significant variation in the duration and lu- 
minosity of the ULIRG phase depending on the local en- 
vironment. Pair mergers would be ultraluminous in the IR 
for a shorter period than multiple mergers and have a lower 
total IR luminosity on average, as there can be more star- 
bursts in a multiple merger. It is worth noting that the com- 
bination of violent interaction coupled with large quantities 
of gas driven to the centres of mergers would undoubtedly 
increase the chances of forming an AGN and hence an opti- 
cal QSO. We argue however that the S88 picture represents 
only a subset of the total ULIRG population described by 
our model. 

To illustrate our alternative model, we describe below 
three points along this continuum of IR power source as a 
function of merger stage. Here we only consider pair mergers, 
as currently no numerical simulations have been performed 
of multiple galaxy mergers (as may be found in compact 
galaxy groups rather than in the field) with varying bulge 
components. 

The first point would be a merger between two bulge- 
less disk/halo galaxies. Starburst activity would be triggered 
early in the interaction while the galaxies are still separate, 
largely exhausting the gas supplies in both galaxies and lead- 
ing to a physical merger that was not ultraluminous in the 
IR due to a starburst. The early starburst in the progenitors 
would make accretion commence early in the interaction, 
possibly leading to an AGN before the galaxies have started 
to physically merge. The merger would be IR luminous for 
the first 10* years or so, followed by a quiescent period be- 
fore (possibly) ending in an IR luminous AGN phase when 
the galaxies merge. The second point would be a merger be- 
tween a disk/halo galaxy and a disk/bulge/halo galaxy. A 
starburst would be triggered in the bulgeless galaxy early in 
the interaction, and will dissipate before the galaxies start 
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to physically merge. After about 5 x 10* years, as the galax- 
ies began to physically merge, a starburst would be trig- 
gered in the disk/bulge/halo galaxy which had been sta- 
bilised against a starburst in the earlier stages of the merger. 
AGN activity could be triggered both early and late in a 
merger between two such galaxies. The third point would 
be a merger between two disk/bulge/halo galaxies. There 
would be no early starburst and hence no ULIRG activ- 
ity whilst the galaxies are still separate. Once the galax- 
ies start to coalesce a strong starburst would be triggered. 
Accretion can only start when the galaxies have started to 
physically merge. IR luminous AGN activity would probably 
commence shortly before the end of the starburst phase. 

This alternative model resolves many of the problems 
with the S88 picture highlighted by previous observations. 
Indeed several observations which directly contradict the 
S88 picture are in excellent agreement with our model. For 
example the study showing that the fraction of ULIRGs 
powered by AGN does not increase with advancing merger 
stage is in strong agreement with our model as our model 
predicts that no such trend should exist. Under our alterna- 
tive model, starburst and AGN activity could be observed at 
any point during a merger, depending on the morphologies 
and local en vironment of the merger p r ogenitors. As previ- 
ousl y no ted (Mihos & Hornquist 1994; Mihos & Hernquist 
199 q) tp is classification scheme would be largely indepen- 
dent of orbital dynamics. 

Although it is not possible with imaging data alone 
to assign our sample to points along this continuum, we 
can give examples of typical morphologies at the points 
described above. A typical 'point 1' object may resemble 
IRAS 02587-6336. Here the galaxies are still widely sepa- 
rated and hence early in the interaction. There is no QSO 
activity but the object is ultraluminous in the IR. A 'point 
2' object might resemble IRAS 20037-1547. A QSO has been 
formed in the bulgeless galaxy which formed an AGN early 
in the interaction due to the early starburst. The source 
IRAS 18580-1-6527 is a good example of a 'point 3' ULIRG, 
where violent star formation has only commenced once the 
galaxies have started to merge. 

An evolutionary track that lasts 10* — 10^ years (such 
as the S88 model or our alternative model) is also testable 
using other quantities that vary on similar timescales, or 
properties that should remain invariant. Observational tests 
include galaxy morphologies, stellar population mixes and 
Mpc scale environments. In particular, the QSOs in our sam- 
ple showed negligible evidence for preferred timing late in 
the interaction though the number of QSOs in our sam- 
ple is small. Increasing the number of ULIRG QSOs se- 
lected purely on the basis of their FIR emission and ob- 
served with very high resolution imaging would provide bet- 
ter constraints on QSO timing in ULIRGs. Possible theoret- 
ical tests include more comprehensive N-body simulations of 
both pair and multiple galaxy mergers, incorporating galaxy 
morphologies, starburst triggering and QSO formation. 



6 CONCLUSIONS 

We have presented a sample of 23 ULIRGs, taken with 
WFPC2 using the F606W band filter on HST, and com- 



pared this to archival HST F814W band data for the same 
objects. Our conclusions are: 

1) There are no isolated normal spirals or ellipticals in 
the sample; the sources are either galaxy pairs, merging or 
interacting systems, or QSOs. The fraction of observed inter- 
acting sources in the sample is 87%. As a whole the sample is 
entirely consistent with ULIRGs being strongly linked with 
both galaxy interactions and (to a smaller extent) with QSO 
activity. In three cases, the HST images showed clear signs 
of interaction which had only been ambiguously detected 
from the ground. Although the majority of the sample are 
physically merging systems a significant fraction (~ 40%) 
exist in binary systems or QSOs. 

2) Stellar population synthesis models indicate that in 
most cases the optical emission is dominated by stellar popu- 
lations with ages of several Gyr or less. Although uncertain- 
ties due to unconstrained mass normalizations and dered- 
dening vectors are significant, it seems likely that any opti- 
cal starburst triggered by interactions does not contribute 
significantly to the optical emission. The average absolute 
magnitude for those sources that do not display a QSO is 
0.3 magnitudes less bright than the field galaxy magnitude 
in the Cousins / band. 

3) Most of the merging systems in the sample show a 
number of compact bright 'knots' whose colour differs sub- 
stantially from the surrounding regions. Colour maps for 
nine of the objects show a non-uniform colour structure in 
all but one case. Observed features include blue regions lo- 
cated towards the centres of merging systems, and compact 
red regions. Although difficult to interpret with only imag- 
ing data the blue regions are likely to be regions of enhanced 
star formation and some of the red regions dust shrouded 
starbursts or AGN. 

4) There is no apparent correlation between the IR lu- 
minosity and either the optical emission or the merger phase 
for non-QSO ULIRGs. Those ULIRGs that contain a QSO 
are on average no more symmetrical than those that do not, 
implying that AGN activity does not tend to reside in more 
advanced mergers. The relatively small numbers of ULIRGs 
with a QSO suggest that this phase is short, and although 
this may be biased by the sample's pre-selected luminos- 
ity distribution, the higher IR luminosity QSOs are if any- 
thing over-represented compared to their relative comoving 
space density. The relatively small number of physically sep- 
arate objects also suggests that this phase is short and that 
ULIRGs are dominated by physically merging systems. The 
QSO host galaxies were found to be either interacting sys- 
tems or ellipticals. Comparisons between our morphologies 
and previous merger models suggest that ULIRGs can be 
both pair mergers between field galaxies and multiple merg- 
ers between galaxies in compact groups. 

5) Based on the morphologies of this sample, previ- 
ous IR spectroscopy, and the results of N-body simulations 
we propose a new model for ULIRG evolution. We pro- 
pose that ULIRGs are not a simple transition stage between 
galaxy mergers and QSOs and there exists no sequence from 
'cool' ULIRG to 'warm' ULIRG to QSO. Under our model, 
ULIRGs are an extremely diverse class of object where the 
time evolution of starburst and AGN activity is driven solely 
by the morphologies of the merger progenitors and the lo- 
cal environment. This model suggests that the diversity in 
ULIRG morphologies and the lack of preferred timing of 
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AGN activity in more advanced mergers is due to a contin- 
uum of merger progenitor morphologies with varying bulge 
components. The local environment (field or compact group) 
affects the number of merger progenitors and hence the fre- 
quency and luminosity of starburst events. Although the 
combination of violent merging and large central gas con- 
centrations would lead to an increased likelihood of AGN 
activity and as such some ULIRGs would represent such a 
transition stage, this is only a subset of the ULIRG popula- 
tion as a whole, and there exists no automatic evolution to 
an optical QSO stage after the ultraluminous IR stage. 
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